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Abstract 

Lieb and Schupp have obtained, using certain version of "spin-reflection positivity" method, a 
number of ground-state properties for frustrated Heisenberg models. One group of these results is 
related to singlet nature of ground state and it needs an assumption of reflection symmetry present 
in the system. In this paper, it is shown that the result holds also for other symmetries (inversion 
etc.). The second Lieb-Schupp result is relation between ground-state energies of certain systems. 
In the paper, this relation is applied to multidimensional models on various lattices. 

PACS numbers: 75.10.Jm; 75.50.Ee; 05.50+q 
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FIG. 1: Example of frustrated spin system. 



I. INTRODUCTION 

Geometric frustration takes place, when no arrangement of spins on the lattice is possible 
in such a way that all interactions minimize their energy (see FIG. QJ. The canonical 
example is an antiferromagnetic Ising model on triangular lattice Another systems 
where the geometric frustration is particularly strong, are antiferromagnetic systems on 
kagome, pyrochlore, or (in d=3) fee cubic lattices. 

Frustrated systems are both very interesting and hard to analyse and understand, both 
in classical version and especially in the quantum case. The source of these difficulties traces 
back to the large ground-state degeneracy in the classical version. Such systems are very 
sensitive to perturbations. A consequence is possibility of very complicated phase diagram 
at finite temperatures. This situation takes place, for instance, in an ANNNI model jsj 
(infinite number of phases, devil's staircase, etc). Besides of numerous efforts and important 
results, 0|, (for reviews, see 0], j^j], 0]), full treatment of such systems is not worked 
out so far. 

The situation for quantum frustrated antiferromagnets is even less clear. It is generally 
suspected that the ground state emerging as a linear combination of many classical config- 
urations is a featureless, "spin liquid" state, i.e. state without long-range ordering, where 
correlation functions fall-off exponentially. However, one cannot exclude another scenario: 
"order by disorder" - exotic orderings absent in a classical version of these models. Such 
scenarios are moreover sensitive to the underlying lattice structure. (For a review, see for 



instance p|). To my best knowledge, no general definite conclusions have been obtained 
here. 

Most informations obtained comes from analysis using approximate methods, such as 
guessing of ground-state wave functions, numerical diagonalization of small systems, spin 
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waves, semiclassical (i.e. large-S) approximations. However, it is difficult to estimate how 
reliable are these approximations. For this reason, exact results are very desirable. However, 
they are very exceptional. 

n such a situation, exact results obtained by Lieb and Schupp , (summarized also 



m 



for fully frustrate systems are of first importance. They are interesting both as 



themselves and moreover, they can serve as a test of validity of approximate methods. Re- 
sults obtained in can be divided into two groups. The first one concern ground-state 
properties for such systems and it can be summarized as follows: ground states are singlets; 
zero-field magnetization is zero; susceptibility is bounded by certain constant. One of as- 
sumptions of these theorems is reflection symmetry. It was apparent that this method could 
work also for another kind of symmetry, not necessarily reflection one. Such a generalization 
turned out to be possible, and this is the first result of the paper. 

The second collection of results concern comparison of ground state energies for spin 
rings. It turned out that similar results hold for more complicated systems (multidimensional 
lattices). This is second result of the paper. 

The outlook of the paper is as follows. In Sec. [Hj the basic ingredients of Lieb and Schupp 
technique is presented. Some systems not mentioned in , [n]] , for which they results hold, 
are listed. Then modifications of their technique to systems exhibiting symmetry other 
than reflection one (rotation or inversion) are described. In Sec JIIIl Lieb and Schupp results 
concerning comparison of ground-state energies for spin rings with the use of certain matrix 
inequality are described. Then it is discussed how this technique can be applied to other 
systems - for instance g.s. energies of systems on different lattices. Section Sec. IIVI contains 
short summary and conclusions. The Appendix contains proof of Lieb-Schupp inequality, 
which is fundamental for Sec. IIHI 



II. LIEB-SCHUPP APPROACH AND GENERALIZATIONS FOR SPIN SYS- 
TEMS WITH SYMMETRIES 

A. Reflection-symmetric systems 

1. Assumptions 

We make the following assumptions concerning systems under consideration. 
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FIG. 2: System exhibiting reflection symmetry. 

1. We consider Heisenberg models for arbitrary spin. (Below we consider mainly the 
s = 1/2 case, but generalization to other spin values is straightforward). 

2. The system is invariant with respect to the reflection about the O axis (d—1 plane for 
d- dimensional system). The system consist of two parts ("Left" and "Right" ones), 
which are interchanged under reflection (see FIG. EJ). The Hamiltonian is a sum of 
three parts: Hl, Hr and He] Hl acts only on the left part (i.e. Hl = h®l), whereas 
Hr acts on the right one (i.e. Hr = 1 ® h). He is the Hamiltonian for "crossing 
bonds" (i.e. for bonds which intersect the symmetry line). It has the form: 

H c = ^2 ^' Si ' Si ' W 

i€L, i'eR 

where is a diagonal matrix with non-negative diagonal elements. Sj can have 
various nature; it can be a single spin or, more generally, it can have the form: Sj = 
^2 aeL ji;a s a, where are real coefficients. The objects Sj can be different for different 
indices i. Under reflection operation, Hl is transformed to Hr and vice versa, whereas 
He transforms into itself. 

3. Hl and Hr are (almost) arbitrary: the only limitations are: both commute with the 
total spin operator, and 

4. moreover, operators are real in the S 3 basis. 
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2. Ground state properties 



Every state of the system can be expressed as a linear combination: 

*P = J2 C ^a®(^U, (2) 

where states ip^ form a real orthonormal base of S 3 eigenstates for the left subsystem, 
whereas (ipp) TO t are corresponding rotated states for the right subsystem. (The rotation - 
meant as an operation in the spin space - acting for the state: \s, m) produces the state 
(_l)«-»»| Sj _ m )). 

The eigenvalue problem: Hip = Eip can be written as a matrix equation for c matrix: 

h LC + C {h R ) T f c (f) T = ^ ( 3 ) 

8=1 y 

where {hL) a p and (hn) Q /3 are real, symmetric matrices for corresponding terms of the Hamil- 
tonian, whereas ty are real matrices, defined for spin operator s y = ^2 aeL jy,a s a (i n the y-th 
"bond") by: = (ipa\s { y' 3) \ipp) and t®^ = i(^\s y 2) \^). Notice the total "minus" sign 
before third term in Q; it is so because replacement of by s y > gives sign change for com- 
ponents 1 and 3, whereas the i factor in definition of t^ gives minus sign for component 
2. 

The energy expectation value, expressed in term of the c matrix, is 

(ip\H\i;) = Tr {cc ] h L ) + Tr (cW) - ^Tr cHf c{tf)K (4) 

y,i 

Lieb and Schupp have established a number of results concerning ground states of reflection- 
symmetric systems. There are: 

1. The c matrix is hermitian, i.e. without loss of generality we can express eigenstates of 
H by hermitian matrix c = c\ (It is proved using the left-right symmetry and a form 
P|) for energy). 

2. Let the c matrix corresponds to the ground state; then, also the |c| matrix corresponds 
to the ground state. Here |c| = v^c?; we take the unique positive square root. (It is 
proved by writing down (jlj) in the diagonal basis). 

3. This implies that we can choose the base for ground states as a collection of functions 
built up from positive semi-definite (p.s.d.) coefficient matrices. 
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4. The overlap of the ground state with canonical spin-zero state (given by unit matrix 
in a basis of eigenstates in either subsystem) is non-zero. 

5. It implies that there exist ground state with spin zero, and moreover, that coefficient 
matrix for this state is positive semidefmite. It follows from the inequality: > E 
(Eb is an energy of the system in the magnetic field b), which in turn is proved using 
the matrix inequality (JHJ). 

6. We have "ice rule" for frustrated units, i.e. for collection of spins appearing in the 
expression Sj • s$/. More precisely, if we denote: Sj = J2aeL^a s a, then the expectation 
of the third component of sites in each crossing bond vanishes for any ground state 

Under additional assumptions, these results can be made stronger. For instance, if the 
system has periodic boundary condition in at least one direction, then all ground states are 
singlets. Moreover, the magnetic susceptibility is bounded from above, both in the ground 
state and in positive temperatures. 

Remark. Lieb and Schupp in and Q have developed general method and (almost) 
haven't given examples of systems where their method could apply; the only example dis- 
cussed is the checkerboard pyrochlore lattice. Some further examples of physical interest, 
where their results are applicable, are: 

• J\ — J 2 2d Heisenberg model: 

h = Ji Sj ■ Sj + j 2 ^2 s i 



considered by numerous authors (see, for instance, |12| and references therein). This 
model is very interesting, as when Jij J\ is varied, it exhibits a transition from ordered, 
antiferromagnetic state to the 'glassy', non-magnetic, spin-liquid-like one. Lieb and 
Schupp results apply for J 2 < \ J\. Most of approximate methods assume that ground 
state is singlet. Lieb and Schupp' results can be used to justify this assumption. 
Moreover, it can serve also as a test of these methods by supplying rigorous upper 
bounds for susceptibility. 
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• Axial Next-Nearest-Neighbour Heisenberg (ANNNH) model [5j]. This is Heisenberg 
model with two coupling constants; we have isotropic coupling J\ between nearest 
neighbours, and moreover, there is coupling J 2 between second neighbours along one 
of the axes (say, z axis). Such models have been used to describe helical and incom- 
mensurate configurations, and Lifshitz points in magnets (j^ and references therein). 
Lieb-Schupp results can be applied when both constants are antiferromagnetic, or 
when J\ is AF and J 2 is F (in this case, one should take the reflection plane to be 
parallel to the line formed by J 2 couplings). 



B. Generalization to other symmetries 

In the course of proofs in the previous Subsection, geometric properties of systems ex- 
hibiting reflection symmetry were not employed. Only assumption which was used was that 
Hl transforms into Hr and vice versa; particular nature of this transformation was not es- 
sential. It suggests that more general symmetry operations than reflection are allowed. It 
is the case; more precise formulation is as follows. 

We make assumptions identical as 1., 3. and 4. in Subsec. Ill A II ; the assumption 2. is 
changed into the following one. 

2'. The system again can be divided by two identical parts "L" and "R", each of them 
is described by the Hamiltonians Hl — h (g) 1 and Hl — 1 <8) h, respectively. He is the 
Hamiltonian for "bonds" between spins in L and R subsystems. It has the form: 

ieL, i'eR 

where the symmetric 7^/ matrix is positive definite. One can view on this property as a 
demand that He has to be "globally antiferromagnetic", i.e. some coupling constants 7^/ 
can be negative (ferromagnetic), however, the whole matrix 7 has to be positive definite. 
The whole system is invariant with respect to the some symmetry operation T, such that 
T 2 = Id; T transforms the "L" subsystem into "R" one and vice versa. The following 
operations: reflection, C 2 rotation (see FIG. EI), or inversion can serve as examples of such 
operation. Under action of T, Hl transforms into Hr, Hr into Hl and He into itself. 

The system described above can be easily transformed into form equivalent to the one 
considered in the previous Subsection. Let us notice that H c is bilinear in Sj, s^. Then, the 
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FIG. 3: System exhibiting inversion symmetry. 

form: He = XligL i'eR 7v' s * ' s *' can ^ e diagonalized by suitable linear transformation in the 
spin variables on "L" and the same transformation on "R". After diagonalization, He takes 
the form 

He = Ji,rcri ■ or 

ieLj'eR 

where the matrix Jj p is diagonal and has only non-negative elements. Then, we have 
system in the form analogous as ( and one can repeat all considerations from previous 
Subsection, obtaining analogous results. One should only remember that the "ice rule" 
concern spin variables after diagonalizing transformation, i.e. we have 

(V>o|c/ + cr/'lV'o) = 0. 

III. TRACE INEQUALITY AND COMPARISON OF GROUND-STATE ENER- 
GIES FOR VARIOUS SYSTEMS 

Lieb and Schupp have proved beautiful inequality for traces of matrices. It is crucial in 
further considerations. 

Theorem (jiol]. [ill]). For any square matrices c,M,N it is true that 

\Tr J McN*\ < ^(Tr c L Mc L M ] + Tr c R Nc R N*), (5) 

where cr = v^c+c, cl = Vcrf are unique (positive) square roots from positive definite 
matrices c^c and cc^. 
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For the convenience of reader, the proof is supplied in the Appendix. Now let us notice 
that the inequality (0) still remains valid if matrices c, M, N are rectangular ones: If c is 
m x n matrix, then M is m x m matrix, N is n x n matrix, u is partial isometry, whereas 
cl and c R are positive matrices of dimension m x m and n x n, respectively. 

Now, let us consider the system without reflection symmetry, which however still consists 
of the "left" subsystem L described by Hl and the "right" subsystem R described by H R . 
The full system will be denoted by L — R, and its Hamiltonian Hl-r contains also the 
"crossing bond" term H c , describing "interaction" of these two parts: H l _r = H L + H R + 
He- Every state of the system can be written in the form (J2J). We change slightly the 
notation and write an analogon of (J2J) as: 

a,P' 

to distinguish between indices referring to the L part (unprimed) and R (primed ones). 
Remember that the c matrix is now, in general, the rectangular one. 
Consider first the situation, where He contains only one bond: 

H c = s • S, 

where s belongs to L and S belongs to R. One can write expression for the mean value of 
energy in the manner analogous as Eq. (J1J: 

3 

E L _ R = Ticc ] h L + Tr c t ch R - ^ Tr [cH^c(T^] , (7) 

where: Jil, Hr are matrices of corresponding hamiltonians. t refers to the left subsys- 
tem, whereas T to the right one; t, T matrices are given by: t^^ = (V'ak^'^lV'a)? 
C = *<^ (2) |^>, and similarly for T: = <^|^- 3 )|^>, T% = i{^\S^f f ). 

Inequality (0) applied to the last term of (J7J) gives: 

3 1 3 

-^TrcHMctTMy > -L ^(Tr c L t^c L (t^ + Tr c R T^c R (T^) 
fi=i fi=i 

Then, we can write: 

1 3 

E l .r > Tr cJh L + Tr JchR - - J](Tr c L t^c L (&^ + Trc R T^c R (T^) 

11=1 
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= Tr c L c L h L + Tr c R c R h R - 1 ^(Tr c L t^c L {&^ + TV c R T^c R (T^) 

1 3 

= -(Tr c L c[/i L + Tr c\c L h L - J^Tr c^c L (^)t) (8) 

+-(Trc fl 4/i fl + Trcic^ - £ Tr c R T^c R (T^). (9) 

How can we interpret two last expressions? They resemble very much ground-state energies 
for the following systems: the first one consists of two copies of L subsystem with " interac- 
tion" Hamiltonian He = s-s', where s belongs to L, whereas s' to its twin copy (let's denote it 
as the L — L system) . The second one consists of two copies of R subsystem with He = S ■ S' 
(it will be denoted as R — R system). More precisely, they are energies of trial functions, 
built up from matrices and c R , respectively. From variational principle, they are greater 
than true ground-state energies, so we have the general inequality for ground-state energies 
for three systems L — R,L — L,R — R: 

2E L _ R > -Ex-L + E R . R . (10) 

Considerations above concerned the situation, where "interaction" part was the only "bond" 
s ■ S. Generalization to multi-bond case is immediate. Let's have general "interaction" 
hamiltonian: 

H c = ^2^2 Ji i' Si ' S J' ; 

i£L j'eR 

where are positive numbers. Then, the inequality ()10|) is still true for 

systems L — L and R — R with "interaction" hamiltonians: He = XlieL {^2j'eR "A?') s * ' s 'i 
for the L — L system, and H c = J2j'eR (SieL -A?') ' ^'j 1 forR — R system. 

Lieb and Schupp have applied inequality (J 10)) to spin rings, getting the following relation 
between ground-state energies Ek for rings with k spins: 

2-En+m > E 2n + Eim 

However, it seems that it can be used in much more general situations, if we consider the 
multi-bond case. 
Examples. 
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FIG. 4: Illustration of inequality (jllj) in the case of spin chains for n = 2, m = 3. 
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FIG. 5: Illustration of inequality for 2d lattices; to = 2, n = 3. 
1. For spin chains (FigEJ), we have analogous inequality: 



(11) 



2. Rather obvious is generalization of this result for systems defined on subsets of Z d 
for d > 1 (rectangles, parallelepipeds etc) (d = 2 case is illustrated on FIG. EJ); as a 
result, we again have inequality where m, n are lengths of systems in direction 
perpendicular to reflection line (plane for d = 3). 

3. Consider system defined on a rectangular subset of square lattice (see FIG. EJ). Divide 
it into L and R part in nonsymmetric manner (in this example, this division is realized 
by the "snaky" line). Then, we obtain the inequality (fTUJ) for systems as pictured on 
FIG. El 

4. Consider "zig-zag ladder" (see FIG.[7|). Divide it into parts L and R by the line going 
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FIG. 6: Division of the system defined on a rectangular subset of square lattice into L and R 
subsets by "snaky" line. Single lines denote coupling constants equal to J, double lines correspond 
to 2J couplings, and triple ones - to 3J. The same convention is used on two following pictures. 



FIG. 7: "Zig-zag" and ordinary ladders. 

through the middle of the ladder and parallel to it. Then, the inequality (jlOj) gives 
relation between ground-state energies of the "zig-zag" ladder and ordinary one, with 
suitable relation between coupling constants (see FIG.[7J). 



5. The same construction can be applied to the "pyrochlore" ladder (see FIG. | 
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FIG. 8: "Pyrochlore" and ordinary ladders. 

6. Ground-state energy relations for systems described in Example 4 can be repeatedly 
used to obtain the inequality between systems defined on triangular and rectangu- 
lar lattices. Namely, consider the system defined on isotropic triangular lattice with 
coupling constant J and with periodic boundary condition in both ( "horizontal" and 
"vertical") directions. Assume, moreover, that in the horizontal direction we have 2 k 
sites and M in the horizontal one (k,M are arbitrary). Let division of the system 
onto L and R parts will be realized by "vertical" plane in such a way that both parts 
contain equal number of sites. Then, by subsequent use of inequality (fTUj). we obtain 
relation between ground state energies of systems on triangular lattice E tr and on 
rectangular one -Erect: 

Etr > -^rect 

where system on rectangular lattice has coupling constants equal J in the vertical 
direction and 2 J in the horizontal one. 

7. Analogous construction can be applied to obtain the ground-state energy inequality 
between systems defined on lattices: isotropic pyrochlore one (E pyro ) and rectangular 
one (E r ect): 

Epyro ^ E rec t 

where system on pyrochlore lattice has coupling constants all equal to J, whereas those 
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on rectangular lattice has couplings J in the vertical direction and 3 J in the horizontal 
one. 



IV. SUMMARY 

In the course of paper, some extension and applications of Lieb-Schupp approach have 
been given. The first result concerns relaxing of demand of reflection symmetry, present in 
their papers. It turned out that also other kinds of symmetry are allowed (inversion, C2 
rotation). The second result is an application of Lieb-Schupp inequality, relating ground- 
state energies of various systems, to numerous new (to my best knowledge) situations. The 
relation between ground-state energies of models on triangular and rectangular lattices can 
serve as an example. 

In the paper jll| Schupp wrote: "There is no doubt that the scheme can be further 
generalized". I consider my paper as a step in this direction, but of course possibilities of 
the Lieb-Schupp scheme seem to be far from exhaustion. 

APPENDIX A: APPENDIX 

Here we supply proof of the inequality 

Proof. Using polar decomposition theorem, we can express the u matrix as: c = ucr, 
where the u matrix is unitary. We have: {ucrv)) 2 = uc^cu^ = {ucr)(crv)) = (ucr)(ucr)^ = 
cc^ = c\, then, because the (positive) square root is unique, we have: ucrv) = cl- 
Analogously, for arbitrary analytic function / defined on positive real numbers, we have: 
uf(cR)u^ = /(cx,). In particular, u^/cr = ^/clu, which implies: c = a/clMa/cr. Now, let: 
P = v) \fclM y/c£u and Q := ^/crN^ ^/cr] then, we have: 

|Tr c ] McN^\ = \TiPQ\ < - (Tr PP 1 " + Tr QQ f ) 

2 

= i(Tr c L Mc L M^ + TtcrNcrN^), (Al) 
where we have used Schwarz inequality for matrices: 

|TrPQ| = lY^PM < ^E(l^| 2 + l^l 2 ) = ^(TrPPUTrQQt). 
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